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The intrinsic Josephson effect in the high-T c superconductors is studied using the layered two- 
dimensional t-J model. The d.c. Josephson current which flows perpendicular to the t-J planes 
is obtained within the mean-field approximation and the Gutzwiller approximation. We find 
that the Josephson current has its maximum near the optimum doping region as a function of 
the doping rate. 
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Since the discovery of high-T c superconductors, there 
have been revealed many novel superconducting prop- 
erties which can not be expected for conventional su- 
perconductors. One of these is the d-wave symme- 
try of the pair potential giving rise to novel properties 
of phase coherence which has been confirmed by both 
experimental and n thjeco;ct^^f^udics in tunneling and 
Josephson effects.trBB'QQ&OBlj'tiJ The strong two di- 
mensionality and the anomalous transport along c-axis 
direction have brought about another important fea- 
ture of the existence of the .intrinsic Josephson effects 
in Bi 2 Sr 2 CaCu20 8 (BSCCO)My) It was reported that 
a small BSCCO single crystal behaved like a series of 
Josephson junctions. It is natural to consider a BSCCO 
superconductor as a stack of superconducting sheets con- 
sisting of CuO bilayers separated by BiO and SrO lay- 
ers acting as insulators. After the discovery of intrin- 
sic Josephson effect by KleinerJiiEj) there are several 
evidences which support the existence .of Jnsfiphspn cou- 
pling between adjacent CuO bilayers. Ej'oli§li3'[LZl ) Up to 
now there are several theoretical works about this prob- 
lem from the view point nf the layered structure.ES How- 
ever, existing theoriesEj'EJElP do not treat the two im- 
portant facts, i.e., short coherence length and the strong 
correlation on the equal footing. Taking account of the 
short coherence length effect, it is suitable to express 
BSCCO by the lattice model. One of the authors calcu- 
lated the intrinsic Josephson current in layered attractive 
Hubbard model within mean field approximations where 
the s-wave conventional paring is considered£3' Later, 
Schmitt et alE$ calculated the Josephson current in the 
same model based on the Quantum Monte Carlo Calcula- 
tion and established the existence of intrinsic Josephson 
current from the numerical calculations with much more 
accuracy. However, in these two papers, strong correla- 
tion effect and the resulting <i-wave symmetry were not 
considered. To overcome this problem, in this paper, we 
study the intrinsic Josephson current using the layered 
two-dimensional (2D) t-J model, which includes the im- 
portant featureSj^bove mentioned. 

The t-J modeES is one of the promising models which 



explain the low-energy excitations in high-T c supercon- 
ductors. Although the analytic solution of this model 
have not yet been obtained, phase diagrams as a func- 
tion of doping rate S and a superexchange interaction—/ 
are numerically studied. .at *X — for one-dimensionO 
and two-dimensionB3'E3E§Ej) Especially, in the 2D t-J 
model, the obtained phase diagram as a function of dop.- 
ing is consistent with actual high-T c superconductors .E2P 

In this letter, we study the intrinsic Josephson effect 
in the layered two-dimensional t-J model. The intrinsic 
Josephson current which flows perpendicular to the t-J 
planes is calculated and the current-phase relations are 
obtained for several doping rates. The Josephson current 
obtained as a function of the doping concentration S has 
a maximum near the optimum doping region. 

The Hamiltonian of the m-th t-J plane is written as 
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where the Hilbert space is defined on the subspace with- 
out double occupancy, i,j refer to planar sites on a 
square lattice, and c™^ , c™ represent fcrmion operators 
within the m-th layer. The spin and the number opera- 
tors in this layer are defined as 



ST = 



mf / - \ m r, 



ol,(3 



A, 



respectively. There seems to be an agreement in the 
conclusion that the c-axis conductivity in the normal 
state generally cannot be explained as a coherent in- 
terlayer transport. We assume, however, that the co- 
herency among the layers is restored in the supercon- 
ducting state. Thus we allow for the electron hopping 
between the m-th and the m + 1-th layers, 



H? = 



-t (c m+1 K n 
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Now the total Hamiltonian of the layered t-J model is 
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written as 



(3) 



mentum space as 

Weff = - E TkC tc 



N s e 



where fi is the chemical potential of the whole system. 

For this Hamiltonian, we consider Gutzwillcr-type 
variational wave functions P<s|$) with Pq = II m) j(l — 
n™n™) being the Gutzwiller projection operator which 
excludes the double occupancy, and |$) being a one-body 
mean-field wave function. Due to the Gutzwiller projec- 
tion, it is usually difficult to calculate the variational en- 
ergy. Thus, we use a Gutzwiller approximation in which 
the effect of the projection is represented by statistical 
weight factors. For the two-dimensional t-J model, it 
is known that this approximation and the variational 
MonteJIJarlo simulation give very similar variational en- 
ergies. c3 

In the Gutzwiller approximation, the expectation val- 
ues of the terms in j are estimated as 
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where (• ■ •) and (■ ■ -}o represent the expectation values 
in terms of Pg\&) and |3>), respectively The renormal- 
ization factors gt , g s ^j and g^3 are determined by the 
ratios of the probabilities of the corresponding physical 
processes in the states Pg\<&) and $}: 

26 4 

9t = ITs' 9s = (T+*p 9n = h (5) 

where 6 = 1 — n is the doping rate. Since the double 
occupancy is excluded in each layer, we can naturally 
assume that the interlayer hopping term is also estimated 
as 



i±xf 



r)=9t{ 



(G) 



Next, we consider the relation among the supercon- 
ducting phases of the layers. We introduce order param- 
eters in the m-th layer as 
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with t = x and y, i + r denotes the nearest neighbor of 
i in the r direction. Since we consider the d a: 2_ y 2-wave 
symmetry, A™ = —A™. When the Josephson current 
flows along the z-direction, the order parameters satisfy 
the following phase relation: 



A™ = A° exp(imq 



(8) 



Thus, the fermion operator is subject to the following 
condition: 



exp(±2<5 



(9) 



Using the Gutzwiller approximation and eq. (9) , we can 
obtain the effective mean-field Hamiltonian of (3) in mo- 



E [Afe 2D c- fci c fcT + A£ 2D 4 T ct fc J , (10) 



where k = (k x , k y , k z ) and &2D = (k x ,k y ) are the three- 
and two-dimensional wavevectors, respectively, /«2D = 
(k x ,k y ) is the two-dimensional wavevectors, N s is the 
total number of sites, 
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and a Xl a yi a z are the lattice constants. 

This Hamiltonian can be diagonalized using a standard 
Bogoliubov transformation: 



HeS = E { Ek ~ 2gttz Smfc * a * sin 7^ ( A t A k 
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where 



e. = je k + \& k2 »\ 2 , I 



& = --(Tfc + T_ fe ) (12) 



The parameters A°, £, [i are determined by the fol- 
lowing self-consistent equations: 
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Then we obtain the ground state energy per site as the 
function of d>: 



Eo(4,)/N. = ^r^ k -E k ) 
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In the following calculation, we take J/t = Jjy/t = 0.2 
and t z /t = 0.1. 

First, we look at the ^-dependence of Eq((/>)/N s . In 
order to see the variation in the energy clearly, we define 
AE(4>) = (E ((p) - E Q (0))/N s . In Fig.2, AE(<f>) is plot- 
ted against the phase 4> for various values of the doping 
rate <). We can see that AE((f>) is the monotonic in- 
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creasing function of 0, and that AE(<fi) increases with 
8 for arbitrary values of </>. When 8 > 0.25 ( the line 
(e) ), however, E((f>) exceeds the energy of the normal 
state, i.e., the superconducting state becomes unstable 
in the region cf) > O.Sir. Let us discuss this situation a 
little more. From eq.(6), we can interpret g t t z as the 
effective amplitude of the interlayer hopping. Thus, as 8 
increases, the coupling strength between layers increases. 
On the other hand, the gap— amplitude A° linearly de- 
creases with increasing of <5.l1P If the gap amplitude A° 
is smaller than the effective interlayer hopping amplitude 
gtt z , the quasiparticle energy 



(I + | A fc2D | 2 - 2g t t z sin k z a z sin ■ 



(17) 



becomes negative for (f> > cf> c , where <p c is a critical phase. 
Consequently, the superconducting state is no more sta- 
ble as compared to the normal state. 
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Fig. 2. The <f> dependence of the Josephson current j(4>)/t for 
t z /t = 0.1. We take (a) 6 = 0.05, (b) <5 = 0.1, (c) S = 0.15, (d) 
8 = 0.2 and (e) <5 = 0.25. 




Fig 



1. The <p dependence of A(<f) for (1) t z /t = 0.1. We take (a) 
= 0.05, (b) <5 = 0.1, (c) 5 = 0.15, (d) 6 = 0.2 and (e) S = 0.25. 
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Next, we calculate the Josephson current defined as 

m = yj^ow. (is) 

Figure 2 shows the Josephson current-phase relations ob- 
tained for various values of 8. The currents show the si- 
nusoidal (^-dependence as expected and have their max- 
imum values jmax(c>) at (f> = ?. When S < 0.2, j max 
increases with 5. This is because the effective mass of a 
Cooper pair, that is, its localization tendency in a layer 
decreases in rough proportion to 1/|A°|,E3 with |A°| be- 
ing the decreasing function of 6. On the contrary, when 
S > 0.25, j max decreases with increasing of 5 because 
the Josephson current vanishes for <f> > <p c . Thus j max 
reaches its maximum in the region 0.2<£<0.25. In 
order to see this situation in more detail, the doping de- 
pendence of jmax is plotted explicitly in Fig. 3. We can 
see that j ma x has its maximum value at 8 ~ 0.23 in the 
present case. _ 

In the recent experiment Jir it has been reported 



Fig. 3. The maximum Josephson current Jmax(<5)/t plotted as a 
function of & for t z /t = 0.1. 



that the intrinsic Josephson current increases with hole- 
doping (in our notation, 8) when the system is in the 
under-doped region. We note here that, if we do not as- 
sume the doping dependence of interlayer hopping as in 
Eq. (10), i.e., the interlayer hopping amplitude has no 
doping dependence, j(<fi) becomes a decreasing function 
of 8. Since the renormalization of interlayer hopping in 
eq.(10) is a natural assumption to treat the layered t-J 
model based on the Gutzwiller approximation, the dop- 
ing dependence of the intrinsic Josephson current ob- 
served in actual experiments can be used to check the 
validity of the Gutzwiller approximations of the layered 
t-J model. Of course there is only a few experiments 
concerning a doping dependence of the Josephson cur- 
rent so far, so that more extensive comparison between 
the theory and experiments is necessary. 
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In this letter, we have calculated the d.c. Josephson 
current which flows perpendicular to the CuC>2 planes 
in high Tc superconductors using the layered t-J model 
within the mean field theory and Gutzwiller approxima- 
tion. The maximum Josephson current has its maximum 
near the optimum doping rate. This is due to the compe- 
tition between the enhancement of the effective transfer 
energy along the c-axis and the decrease of the mag- 
nitude of the order parameter with the increase of the 
doping rate. We hope such a behavior will be observed 
in experiments near future. 
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